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Abstract 

We present a high-statistics lattice calculation of the kinetic energy — Ai/2mft of the 
heavy quark inside the i?-meson and of the chromo-magnetic term A2, related to the 
B*-B mass splitting, performed in the HQET. Our results have been obtained from a 
numerical simulation based on 600 gauge field configurations generated at (3 = 6.0, on 
a lattice volume 24 3 x 40 and using, for the meson correlators, the results obtained with 
the SW-Clover 0(a) improved lattice action for the light quarks. For the kinetic energy 
we found -Ai = (B\h(iD) 2 h\B) / \2M B ) = -(0.09 ± 0.14) GeV 2 , which is interesting 
for phenomenological applications. We also find A2 = 0.07 ±0.01 GeV 2 , corresponding 
to M|, -M| = 4A 2 = 0.280 ±0.060 GeV 2 , which is about one half of the experimental 
value. The origin of the discrepancy with the experimental number needs to be clarified. 



On leave of absence from Dept. of Physics, University of Southampton, Southampton S017 1BJ, UK. 



1 Introduction 



In a previous study ||] we presented several results for the 5-meson binding and kinetic 
energies, obtained by using the lattice version of the Heavy Quark Effective Theory 
(HQET) 0. The calculation required the non-perturbative subtraction of the power 
divergences present in the matrix elements of the Lagrangian operator hD^h and of the 
kinetic energy operator hD 2 h, following the strategy outlined in ref. ||]. Good results 
were found for the binding energy A, which have been confirmed and improved by the 
higher statistics analysis of ref. On the other hand, given the low statistical sample 
at our disposal in ref. [jj], the determination of Ai was quite poor and resulted only in 
a weak upper bound of |Ai| < 1 GeV 2 . 

In this paper, we present a high-statistics calculation of both Ai and A2, based on a 
sample of 600 configurations at = 6.0, on a lattice volume 24 3 x 40. Using the same 
renormalized kinetic energy operator as in dimensional regularization (e.g. as in the 
MS scheme) Q we find Ai = 0.09 ± 0.14 GeV 2 . Even though the determination of Ai is 
still not very precise, it is however interesting for phenomenological applications, and 
for comparison with other theoretical predictions. Our result for the matrix element of 
the chromo-magnetic operator, A2 ~ 0.07 GeV 2 , obtained with rather large statistics, 
is in excellent agreement with the latest result of the UKQCD collaboration |H], which 
had been performed with a smaller statistical sample, but at a smaller value of the 
lattice spacing 0. A value of A 2 ~ 0.07 GeV 2 corresponds to M§» -M§ ~ 0.280 GeV 2 , 
which is about one half of the experimental number, 0.485 GeV 2 . In our opinion there 
are three possible sources which may contribute to the explanation of the discrepancy. 
The first is the fact that, at one- loop order, the factor necessary to match the chromo- 
magnetic lattice operator of the HQET to its continuum counter part is very large || , 
1 Our kinetic energy operator is made finite by subtracting non-perturbatively the power divergences, 

and then by matching it to the MS operator using perturbation theory. 

The original version of the UKQCD preprint @, had an error of a factor of 2 in the value for A2. The 

corrected value of the mass splitting is M§„ - M% = 0.281 ± 0.015 ± 0.040 GeV 2 . 
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and it is possible that the higher-order terms modify this factor significantly. A second 
source of systematic error is quenching. The final possibility is that the discrepancy is 
due to a physical reason, i.e. that M^* — M\ receives significant contributions from 
higher-order terms in the HQET. It seems to us unlikely that any of the three sources 
alone could explain a discrepancy of about a factor of two between the lattice and 
experimental values of A2. 

The calculation of Ai and A2 proceeds in two different ways. The chromo-magnetic 
operator that enters in the calculation of A2 is logarithmically divergent in the ultra- 
violet cut-off, i.e. the inverse lattice spacing 1/a. Its multiplicative renormalization 
constant can in principle be computed in perturbation theory or with the non- 
perturbative technique proposed in ref. 0]. In the calculation of Ai, instead, one 
has to make a non-perturbative subtraction, related to the presence of a quadratic 
divergence in the matrix elements of hD 2 h ^. The subtracted operator only requires a 
finite multiplicative renormalization constant, which is different from unity because of 
the breaking of the reparametrization invariance of the HQET on the lattice @|— [10|- 

Although it may be convenient to introduce the parameter Ai (or the binding energy 
A), it is not necessary to do so. One can either relate two or more physical quantities 
directly to the required precision in the heavy-quark expansion, or determine physical 
quantities directly from some non-perturbative computation, such as for example from 
the calculation of the bare lattice matrix elements in numerical simulations. Indeed the 
actual values of A and Ai by themselves have no physical meaning, since they depend 
on the definition that one decides to adopt for the corresponding operators, and do not 
give us any relevant physical information. Predictions for physical quantities can be 
obtained only by combining these parameters with the coefficients and matrix elements 
of the dominant operators. This is to be contrasted with a s or A2, which correspond 
directly to measurable physical quantities up to higher power corrections in the inverse 



3 In order to define a finite operator hD 2 h the subtraction of a linearly divergent term is also needed. 
This subtraction, however, does not enter in the calculation of Ai, see below. 
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mass of the heavy quark. The introduction of a consistent definition of Ax, which is 
independent of the regularization and is of 0(Aq CD ), presents us, however, with an 
intuitive estimate of the size of the corrections and can be used for comparison with 
other definitions currently used in heavy-flavour physics. 

The plan of the paper is the following. In sec. || we recall the relevant formulae, 
which define the non-perturbative procedure for renormalizing the operators hD 2 h 
H and the method used to compute the hadronic matrix element of the subtracted 
operator in order to obtain Ai; in the same section we also describe the method used to 
obtain A2; in sees. ||| and |] we present the numerical results for Ai and A2, respectively, 
together with a comparison with other calculations of the same quantities. Finally in 
section || we present our conclusions. 

2 Definition of the renormalized Ai and A2 

In this section we define the prescriptions that we will use to calculate Ai and A2. 
We also recall the procedure used to extract the hadronic matrix elements of the 
renormalized operators from suitable two- and three-point correlation functions. 

The renormalization of hD 2 h can be implemented by imposing suitable non-per- 
turbative renormalization conditions on its quark matrix elements in the Landau 
gauge 

Si- This is not necessary and it would certainly be possible to obtain physical 
predictions using the matrix elements of the bare lattice operator, which diverges as 
powers of 1/a. Indeed, in a theory regulated by a dimensionful cut-off, it is consis- 
tent not to perform the subtractions of the power-divergent terms at all, but to work 
with the bare operators and to compute the coefficient functions (which will therefore 
contain powers of the cut-off) in perturbation theory |3|, [ll]]. In this case, however, 
the matrix elements in the effective theory are divergent in the ultra-violet cut-off 
and depend on the regularization. Therefore they cannot be interpreted as "physical" 
quantities, in contrast with the approach that we have adopted in ref. Q and here, 
and with the use of Ai that is made in the literature. Since Ai enters in the theo- 
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retical predictions for several quantities relevant to the physics of hadrons containing 
the 6-quark (hadron spectroscopy, inclusive decays, form factors in exclusive weak and 
radiative decays, etc.), we find it convenient to introduce a quantity that is finite and 
independent of the regularization |3| . 

In order to remove the power divergences from the kinetic energy operator, we have 
imposed on the relevant operator a renormalization condition that corresponds to the 
"physical" requirement (h(p = 0)\hD 2 h\h(p = 0)) = 0, where hD 2 h is the subtracted 
kinetic energy operator ||]. Since a complete discussion of the main theoretical as- 
pects of the proposed non-perturbative renormalization conditions and matching of 
the operators has been presented in our previous papers on this subject jl], [|, only 
those formulae which are necessary to the understanding of the numerical results are 
reported in the following. 

The subtracted kinetic operator hD 2 h, which is free of power divergences, has the 
form 



h(x) D 2 h(x) = h(x) D'i h(x) - — 

a 

- %h(x)h(x), 



(1 + a 5m) 



(h(x) D^h{x) + 5mh(x)h(x) 



where for the discretized version of -D4 and D 2 we have taken 
Wop = \ ( 5 *P $x,x ~ U^(x - ai)5 x x 



at ' 



and 



D 2 



3 

a- ' 



4 E (U k ap {x)5 xx + a% + U k J(x - ak)5 xx _ ak - 25 aP 8 x>x ) . (3) 
a k=i 

The term in square parenthesis in eq. (|l|) is the subtracted lattice Lagrangian operator 
at lowest order in the HQET and 5m is the suitable mass counter-term needed to 
eliminate the linear divergence present in the unsubtracted Lagrangian []. 



For simplicity we give here the expression corresponding to the unimproved lattice Lagrangian for the 
heavy quark; the more complicated expression necessary in the improved case can be found in ref. S. In 
the following, the numerical results will be given for the improved case only. 
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In eq. (|T|) the constants c\ and C2 are functions of the bare lattice coupling constant 
go{a). They have been computed in one- loop perturbation theory in ref. Notice 
that we have preferred to express hD 2 h in terms of the subtracted Lagrangian operator, 
which explicitly contains the residual mass 5m. In this way we can use the equations 
of motion of the subtracted Lagrangian. This will prove useful below. 

In order to eliminate the quadratic and linear power divergences, a possible non- 
perturbative renormalization condition for hD 2 h is that its subtracted matrix element, 
computed for a quark at rest in the Landau gauge, vanishes: 

(h(p = 0)\hD 2 h\h(p = 0)} = 0. (4) 

This is equivalent to defining the subtraction constants through the relation (in the 
following we will work in lattice units, setting a = 1) 

Pfr^) = c i + c 2*> ( 5 ) 

where 

,.. Ef.oEi, f <sOMI?,0 »?(«') S(S*'|5,o)} 
'* W S E.(g(g.«|8.o» ' <6> 

where y = (y,f). S(x,t\0,Q) is the heavy-quark propagator between the point of 
coordinates (x,t) and (x = 0,t = 0) computed on a single gauge field configuration; 
(...) denotes the average over the gauge field configurations. 

For most applications it is only the constant c% that is required; C2 can also be 
determined directly by eliminating the sum over t' in eq. (||): 

l+l . Zx,y(S(x,t\y, t')D 2 y (t')S(y,t'\0,0)) 
c 2 = Pfrit ,t) = — * — V— (7) 

Ex\S(x,t\o,o)) 

for t' 7^ 0, t. The term proportional to c\ disappears because, by the equations of 
motion of the subtracted Lagrangian, it can only contribute to contact terms. 

The renormalized kinetic energy Ai of the 6-quark inside a l?-meson is then given 

by 

Ax = Z 3l (Aj« - c 2 ) = (l - gx 5f ) (Aj« - c 2 ) , (8) 



where \\™ e = (B\hD 2 h\B) / (2M B ). The term proportional to Xg 2 in eq. (g) is ab- 
sent in continuum formulations of the HQET, and is a manifestation of the lack of 
reparametrization invariance in its lattice version. It has been calculated in ref. ||. 
One can argue that the counter-term C2 defined above is gauge invariant, in spite of 
the fact that it has been derived from quark Green functions in a fixed gauge. As a 
consequence, since the matching constant Zq 2 is gauge invariant, Ai defined in eq. (g) 
does not depend on the gauge. 

^bare can ^ e determined from a computation of two- and three-point i?-meson 
correlation functions in the standard way. Given the quantities 

C B (t) = ^<0|J B (x,i) 4(0,0)|0) and 

X 

c 32 (t',t) = ^<o|j B (x,i) My,0^(OM^04(o,o)|o), (9) 

we have, for sufficiently large values of t' and t — t': 

C B (t) -» Z 2 exp (-(£ - Sm)t) , 

C 32 (t',t) Z 2 X\ am exp {-{£ - Sm)t) , (10) 

where Jg(x, t) (Jb{x, t)) is a source which creates (annihilates) a pseudoscalar B-meson 
state from the vacuum; 8 is the "bare" binding energy and 5m = 5m(5m) a suitable 
mass counter-term ||, ||]. A convenient way to extract Aj" 116 is to compute the ratio 

As usual Aj* are must be evaluated in an interval in which Rg 2 {t',t) is independent of 
the times t' and t, so that the contribution from excited states and contact terms can 
be neglected. 

The computation of the chromo- magnetic matrix element A2, defined as 

A^ 6 = l -(B\h ~ aij G ij h\B)/{2 M B ) , (12) 

proceeds in a similar way. The most important difference is that A2 is free of power 
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divergences so that, at least in principle, there is no need of a non-perturbative sub- 
traction 

\ 2 = Z s . 6 \r\ (13) 

where Z- a is the renormalization constant necessary to remove the logarithmic diver- 
gence present in the bare operator ||. In eq. (|l2|), = 1/2 [7i,7j] is the spin Dirac 
matrix and G %3 is the colour field-strength tensor. In our computation we use the 
clover definition of G lJ introduced in ref. fl2||. To determine A2, we also used a ratio 



of correlation functions as in eq. (11), with the chromo-magnetic operator 1/2 a-G 
inserted in the three-point function. In this case we will denote the ratio (and the 
corresponding three-point function) as R-g(t',t) (C^ g{t' ,t)). 



3 Numerical calculation of Ai 

As explained in the previous section, the renormalization procedure, which defines a 
finite kinetic energy operator, requires the computation of c\ and C2. These subtrac- 
tion constants were obtained using a sample of 600 configurations on a 24 3 x 40 lattice 
at (3 = 6.0. Also A] 3are has been computed on the same set of gauge field configu- 
rations as ci and C2. The relevant two- and three-point meson correlation functions 
were computed using the improved SW-Clover action [|lj| for the light quarks (with 
improved- improved propagators [Q|), in the quenched approximation, at three values 
of the light-quark masses corresponding to K = 0.1425, 0.1432 and 0.1440, where K 
is the Wilson parameter. As in our previous studies of the lattice HQET 113], for the 
heavy meson sources we used the standard axial currents smeared over cubes of size 
Ls = 5, 7 and 9. All the errors have been computed with the jack-knife method by 
decimating 12 configurations at a time. 

The same 600 configurations, quark propagators and smeared sources were used 
for the calculation of both Ai and A2- Thus we will not repeat the discussion of the 
parameters of the numerical calculation when discussing the results for A2. 
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3.1 Determination of the subtraction constant C2 

In order to compute the renormalized Ai, we only need to determine the subtraction 
constant C2, cf. eq. (||). There are several methods to define and extract C2 from the 
heavy-quark Green functions: 

1. One possibility is to fit the time-dependence of pg2(t) at large time distances 
(t — > oo) to a straight line, see eq. (||). This corresponds to imposing the renor- 
malization condition (Q) at p± = 0. Previous studies of the propagator at large 
time distances |], IlL an d the results presented below, suggest that the infrared 
limit of the ratio in eq. (|]) does indeed exist, in spite of the confinement effects 
associated with the heavy-quark propagator. 

2. Similarly, it is possible to define C2 from pg 2 {t' , t), see eq. (^), in the limit t' — > oo 
and t — t' — > oo. 

3. A different possibility consists in defining C2 at fixed t' ^ and t — t' ^ (for 
simplicity we only discuss the case t* = t' = t — t'). We denote the corresponding 
subtraction constant by C2(t*): t* parametrizes the renormalization prescription 
dependence and can be considered as the renormalization point in coordinate 
space. Since the lattice kinetic energy operator has no anomalous dimension, we 
do not expect any dependence of C2(t*) on the renormalization point t* at any 
order of perturbation theory, and this is supported by our numerical results []. 
The use of two- and three-point Green fuctions at small times, t*\QCD 1> to 
define C2 and hence Ai does not require any assumption about the behaviour of 
the heavy-quark propagator at large times. Our numerical results do not show 
any visible dependence (up to contact terms) on t* . 

4. On a lattice with a finite lattice spacing, contact terms are not localized on a 
single point, t' = or t 1 = t, but are smeared over several lattice spacings (see 
for example ref. and fig. ||, to be discussed below). For this reason, besides 

5 This is true when t* is chosen in such a way as to avoid contact terms. 
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Figure 1: The ratio pg 2 (t',t) defined in eq. (Qj, at (3 = 6.0, as a function oft', at several 
values of t, t = 6, 7, 8 and 10. 

Pr>2{t) = J2t'=o P3 2 ^' ''*)' ^ n or der to reduce the effect of contact terms in the 
extraction of C2, we have also linearly fitted the quantity 

t-A 

P% (*) = E P32 ci + c 2 t , (14) 

t'=A 

where the points close to t and to the origin have been eliminated from the sum. 

In fig. [I], we plot pg 2 (t', t), as defined in eq. (0), as a function of t', at fixed t = 6, 7, 8 
and 10. Up to contact terms, we expect pg 2 (t',t) to be a constant in t', at fixed t, 
and also to be independent of t. If the contact terms were entirely due to the mixing 
of the kinetic energy operator with the inverse propagator, as in eq. ([j]), we should 
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c 2 from pg 2 (f,t) and p^ 2 2 (t) 


t 


Interval 


x 2 




6 


3-3 


0.0 


-0.728(2) 


7 


3-4 


0.1 


-0.729(3) 


8 


3-5 


0.1 


-0.729(4) 


9 


3-6 


0.4 


-0.727(6) 


10 


3-7 


0.2 


-0.721(10) 


Interval 


C2 


(6 


: 7)-(10 : 


16) 


-0.727(6) 



Table 1: Values of c 2 obtained as explained in the text. We also give the uncorrelated 

find two spikes, at t' = and t' = t, and a constant value of p^ 2 (t',t) for i' 7^ 0,i. 
The presence of operators of higher dimension, due to discretization errors, introduces 
terms that behave as derivatives of ^-functions (in time), giving rise to the bell-shape 
behaviour of pg 2 (t',t) shown in fig. [I]. Thus in order to obtain c 2 , we have to look for 
a plateau in the central region in t' , at large values of t. ^From the figure, we see that, 
since the contact terms extend over about two lattice spacings in both t' and t — t', 
it is possible to recognize a plateau in t' only for t = 7-10 (at larger values of t the 
statistical errors are too large to draw a firm conclusion). At t = 6 the best that we 
can do is to consider only the point at t 1 = 3. To obtain c 2 we have made a weighted 
average of the values of pg 2 (t',t) over the plateau region, at fixed t. The results are 
given in table [|. In the table, the column Interval denotes the interval in t! over which 
Pnz^'ji) has been averaged. 

If we fit pg2{t) to a straight line on intervals that include small values of t, the 
presence of extended, time-dependent contact terms can also induce a systematic error 
in the determination of c 2 from this quantity. For this reason, in the extraction of c 2 
from a fit to p^ 2 {t) we only used p^ 2 {t) with t > 6. The result is also given in table [j]. 
For this quantity Interval denotes the time interval on which p^ 2 {t) has been fitted 
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and (6 : 7)-(10 : 16) indicates minimum-maximum values of the time on which the fit 
has been made, i.e. the fits were obtained starting with t = 6 or 7 and ending at any 
of the points between 10 and 16. 

The results given in the table show that the values of C2 as obtained from pg 2 (t', t) 
are essentially independent of t (and of t'), as expected from the vanishing of the 
anomalous dimension, and are indistinguishable from the value extracted with a linear 
fit to p*~p 2 (t). Our best estimate for this quantity is then 

c 2 = -0.729 ± 0.005 , (15) 

to be compared with our previous result c 2 = -0.73 ± 0.02 Q. In eq. (||), the central 
value has been taken from Pfivit' = 3—5, i = 8); the error includes the fluctuations of 
pg 2 (t',t) between different points in t', at fixed t, as well as the variation of pg 2 (t',t) 
with t or variations of C2 as derived from a fit to P^ 2 (t) with respect to different time 
intervals. 

3.2 Determination of A^ 6 

The procedure to extract operator matrix elements is standard. It is the same as 
the second method that we have used in the previous subsection to determine c%. At 
fixed t, we study the behaviour of the ratio Rg 2 (t',t) as a function of t', searching for 
a plateau in t'. Here Aj" 116 is defined by the weighted average of the data points in 
the central plateau region, if this exists. We will take as our best determination of 
A] 58 * 6 , the value evaluated in a time interval where the ratio Rg 2 (t\t) appears to be 
independent of both t and t' . In addition, we have to require that the lightest state 
has been isolated. With the smeared sources used in the present case, we know that 
this happens at time distances (t — t') and t' > 4-5. This implies that the total time 
distance t for Rg 2 (t', t) has to be at least 8-10. Moreover, by using (t — t') and if > 4-5, 
we eliminate the contact terms which, on the basis of the discussion in the previous 
subsection, are expected to be present up to distances of order 2-3. 
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Figure 2: 27ie ratto Xf = Rg 2 (f, t) as a function oft' , at t = 8, for smearing sizes L$ - 
and 9. 

In the two- and three-point meson correlation functions, in order to improve the 
isolation of the lightest meson state at short time distances, we have used, as in our 



previous studies [15], the following sources 

Jb{x) = ^h{xi)lols<l(x) 

i 



(16) 
(17) 



where we sum the position of the fields Xi over cubes of sizes L$ = 5, 7 and 9, centred 
on x. 

In order to show the dependence of the results on Lg, we give in fig. Q the results for 
Af D = Rg 2 {t',t) as a function of t' ', at fixed t = 8, obtained using the double-smeared 
sources Jg(x), at the value of the light-quark Wilson parameter K = 0.1432. We 
observe some dependence of Af D on L5. For L$ < 5, this is not surprising, since we 
know, from previous studies of smeared-source correlators, that it is not possible to 
have a good isolation of the lightest state for Lg < 5. Between L$ = 7 and L$ = 9, 
there is a small residual shift in the value of X^ D which will be taken into account 
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Figure 3: The ratio \® D = Rg 2 {t',t) as a function oft', at t = 8 for Ls = 7. The results 
obtained in ref. [1] and those of the present study are shown on the left-hand and right-hand 
side of the figure respectively. 

in the evaluation of the error on Aj ,are . In the following, we will use L$ = 7 with 
double-smeared sources (which was found to be the optimal cubic smearing p5| ) to 
determine the central value of \\ are . 

In fig. |3|, we give, for Ls = 7 and at K = 0.1425, the plateau for \® D as a function of 
t f at t = 8. In the figure, the extracted value of Ai are , obtained by a weighted average of 
Af D in the interval t' = 3-5 (the same as the corresponding interval used to determine 
C2) is also given (solid line). The band limited by the dashed lines corresponds to the 
error of this quantity. For comparison we also report the results obtained in ref. jj]], at 
the same value of (3, on a volume 16 3 x 32, with a statistics of only 36 configurations. 
The significant reduction of the statistical errors with respect to our previous study 
allows for a much better identification of the plateau region, which had been assumed 
to exist, but could not be really observed, in ref. 0]. Morevover, in the present case, 
we have computed the meson correlation functions at three values of the hopping 
parameter K, which allows us to extrapolate in the mass of the light quarks Aj )are to 
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A barc from R 32 (t\t) and Rg 2 (t) for L s = 7 


K 


t = 8, t' = 3-5 


Fit with 6 < t < 16 


0.1425 
0.1432 
0.1440 


-0.74(3) 
-0.74(3) 
-0.73(3) 


-0.70(4) 
-0.70(4) 
-0.69(4) 


0.145431 


-0.71(3) 


-0.67(4) 



Table 2: Values o/A 1 j >are obtained as explained in the text. 



the chiral limit (K c = 0.145431(10)), corresponding to the kinetic energy of the heavy 
quark in the meson. 

Another method to extract Aj )are is to introduce the quantity 

t 



R 32 (t) = J2R32(t , ,t) = a 1 + \ h 1 ™ e t. 



(18) 



t'=o 



and to fit it to a straight line. This method was used in old numerical calculations of 



hadronic matrix elements, such as the kaon i3-parameter [17] and the nuclueon cr-term 
[18], and more recently by the UKQCD collaboration for the calculation of A2 |J, see 
also below. A possible inconvenience of this method is due to the fact that the presence 
of contact terms, and excited states, when t is not large enough can significantly distort 
the value of the slope in time, and hence of the value of the matrix element extracted 
from the correlation function. In this respect, the use of Rg 2 (t',t) is better in that it 
allows us to verify the existence of a plateau, and the independence of the results of t 
and t'. 

We give in table § the results obtained by using the two methods outlined above. 
From this table, we see that the typical statistical error is of the order of 0.03. Taking 
into account the fluctuations of the results between different points in t', at fixed t, and 
in t, of the variations between results obtained with the linear fit to Rg 2 (t) on different 
time intervals, and of the differences between L$ = 7 and L$ = 9, we estimate, for the 
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value extrapolated to the chiral limit: 

A bare _ _Q gg ± Q gg ± q gg ^g) 

In order to appreciate the reduction of the errors, we note that our old result was 
A bare = _ 072 ± 0.14 at K = 0.1425 [|. 

3.3 Determination of the "physical" value of Ai 

We are now in a position to compute the "physical" value of Ai, i.e. the finite matrix 
element of the subtracted operator. In this section the lattice spacing a is written 
explicitly in the formulae. For the physical value of the inverse lattice spacing we use 
a -1 = 2.0 ± 0.2 GeV, where the uncertainty takes into account the spread of values 
obtained in the quenched approximation from different determinations of a (from the 
string tension, from the mass spectrum, from f n , etc.). For the subtraction constant 



C2, we use the result given in eq. (15) 



a 2 c 2 = -0.729 ± 0.005, (20) 



together with 



a 2 A^ 6 = -0.69 ± 0.03 ± 0.03 . (21) 

We thus obtain 

Ai = a- 2 Z 3 2 (a 2 Ai are - a 2 c 2 ) = 0.18 ± 0.22 GeV 2 . (22) 
We have taken Z^ 2 = 1.13 obtained from the expression of the one-loop renormaliza- 



tion constant M 



where 

3 

A = (1 - cos(fcia)) , (24) 



i=l 
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and used the boosted strong coupling constant [19 



uq is a measure of the average link variable for which we take uq = (8i^ c ) _1 with 
K c = 0.145431(10). The value of the renormalization constant in standard perturbation 

theory would be Z^ 2 = 1.07. One can try to reduce the error by taking the difference 

u s 

between Aj )are and the counter-term for each jack-knife individually. In this case we 
obtain 

Ai = 0.09 ± 0.14 GeV 2 , (26) 

which we take to be our best result. 

One could argue that the subtraction is not really necessary, since the effective 
theory on the lattice does not have renormalons. Even though this is indeed true, the 
difficulty in the determination of corrections of order 1 / itiq related to the kinetic energy 
operator would remain the same. The argument goes as follows. The bare kinetic 
energy operator has a very large matrix element a~ 2 x (a 2 Aj >arc ) ~ 2 2 x —0.69 GeV 2 = 
—2.8 GeV 2 , while one expects a correction due to the kinetic energy of the heavy quark 
of the order of the squared Fermi momentum p 2 F ~ Aq CD ~ 0.1-0.6 GeV 2 . Thus the 
huge contribution of the matrix element of the bare operator has to be compensated by 



the corresponding term in the coefficient function of hh [1C]. This requires an extreme 
accuracy in the perturbative calculation of the coefficient function and it remains true 
in the subtracted as well as in the unsubtracted case []. 



In the past Ai has been computed in different models and using QCD sum rules [21 



[23]. Bounds on Ai were also obtained from zero-recoil sum rules [24], which however are 



weakened by higher-order perturbative corrections [25]. The extraction of its value from 
experimental data has also been attempted small values of Ai are suggested 

by the theoretical analysis of ref. |3(| and by the recent studies, ]^0], [HJ. The results 
of the different theoretical estimates are given in table || 



6 Similar problems are encountered also in different approaches to determine Ai , see for example ref. p0[ . 
7 The value attributed to ref. Q m table ||[ was extracted from the results of this paper by Neubert in 
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Reference 



Eletsky and Shuryak [22| 
Ball and Braun ||2 
Bigi et al. J 



Kapustin et al. [25| 



Ligeti and Nir |2q| 



Gremm et al. [ 20 1 



Chernyak [31] 



Method 



QCD sum rules 
QCD sum rules 
ZR sum rules 
ZR sum rules +0(a s ) + O(a^/?o) 
Experiment 
Experiment 
Experiment 
Experiment 



-Ai (GeV 2 ) 



0.18 ±0.06 
0.52 ±0.12 

> 0.36 
No Bound 

< 0.63 if A > 240 MeV 

< 0.10 if A > 500 MeV 

0.19 ±0.10 
0.14 ±0.03 



Table 3: Some of the values of Ai obtained in different theoretical analyses. "Experiment" 
denotes the extraction of X± from the experimental data, for example the charged-lepton 
spectrum distribution in semileptonic B -meson decays. "ZR sum rules" denotes the zero 
recoil sum rules. 

A comparison of the different results and limits given in table [| between them 
and with ours is difficult for two reasons that we briefly discuss. We have seen that 
the subtraction of the quadratic divergence is necessary in order to define a kinetic 
energy operator, the matrix elements of which are finite. The precise definition of 
the renormalized operator, and hence the value of Ai, depends on the renormalization 
procedure and different definitions have been used in the results of table ^. For ex- 
ample in ref. j2(| the definition of Ai is implicitly given by the weight function Wa, 
which depends quadratically on the "ultraviolet cut-off" A, similarly to the quadratic 
dependence of the lattice Ai are on the inverse lattice spacing. The relation between the 
different definitions is in principle computable in perturbation theory, although this 
exercise has not been done to date. The second, related remark (which is general for 
power-suppressed corrections in all effective theories) concerns the precision that can 
be achieved in evaluating, for a given physical quantity, the power corrections at a given 



ref. 



The value for ref. [E0| was taken from the revised version of this paper. 
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order in the inverse quark mass 1/mg. As discussed in ref. [ 10 1 , there are essentially 
two procedures to compute the power corrections. One is based on the determination 
of the parameters of the HQET, such as A and Ai, from the measurement of some ex- 
perimental quantity, as was done for example in ref. |2(|. The parameters found with 
this method can then be used to predict other physical quantities. The other approach 
is that followed in this paper. The parameters of the HQET are defined and computed 
with some non-perturbative technique, the lattice numerical simulations in our case, 
and then used in the calculation of physical quantities, with the same renormalization 
scheme. In both cases, however, it turns out to be rather difficult to achieve a sufficient 
precision in the calculation of the power corrections, because of the truncation of the 
perturbative series for the Wilson coefficients of the 1/rriQ expansion. The truncation 
of the perturbation series limits also the possibility of accurate comparisons between 
values of A and Ai obtained in the different approaches. A detailed discussion of this 
problem, with several examples, can be found in ref. |]l0| (see also Q). The disappear- 
ance of the bound found in ref. |jll| due to higher-order perturbative corrections, which 
are quadratically divergent in the cut-off (~ a s A 2 ) |25( j, is just an explicit example of 
what is expected on the basis of the general arguments of refs. [||, 10|. This holds 



true also in the case of ref. [20|, where the one-loop perturbative corrections to the 
quantities used to obtain A and Ai are very small. For example, when using the results 



of this analysis to predict the mass of the 6-quark in the MS scheme, an intrinsic error 



of about 200 MeV from higher-order perturbative contributions is expected [1C], much 
larger than the correction given by the kinetic energy operator, which is of the order 
of 30 MeV. For the above reasons, although some of the results in table ^ look barely 
compatible with each other and with our result, the only sensible thing to do is to 
compare the predictions of the different approaches for some physical quantity (some 
width, the MS mass, etc.) to the desired accuracy in 1/mg (or 1/wig), and subject 
to the condition that the perturbative corrections are under control. The present dif- 
ferences can be taken as a demonstration that these higher-order corrections are likely 
to be considerable. 
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Although Aj >are for a single hadron is not a physical parameter, the difference of 
the 6-quark kinetic energy of two different hadrons is a well-defined quantity which 
does not require any subtraction. For example, we can study the kinetic energy in a 
l?-meson as a function of the mass of the light quark in the meson, by fitting 

X h r e = A Xl + B Xl M 2 PS , (27) 

where Mp s is the mass of the pseudoscalar meson made by two light quarks with the 
same mass as the light quark in the i?-meson. From our data, by assuming that Mp s 
depends linearly on the value of the quark masses, we get 

Ai(B s ) - Xi(B d ) = -0.09 ± 0.04 GeV 2 , (28) 

where we have also included the factor due to the renormalization constant of the 
lattice bare operator. The result in eq. (|28| ) compares well with the number that can 
be extracted from the following combination of heavy-meson masses: 

1/2 (1/M D - I/Mb) \ rn Q J 

(29) 

where the bar over the meson masses indicates "spin average", e.g. Mp = {Mb + 
3M#*)/4, and rriQ denotes the mass of the c- or 6-quark. This value obtained from 
spectroscopy should be treated with some caution, since its determination involves 
cancellations of much larger quantities, and it is possible that higher-order terms affect 
the result. Similarly one can compute Ai(A^) — Ai(-Bd), by extrapolating both the 
baryonic and mesonic kinetic energies to the chiral limit, and taking the differerence 
(or the other way round). 



4 Numerical calculation of 



The numerical calculation of AJf 1 re proceeds along the same lines as the calculation of 
A^ 6 , so that we will not repeat all the details given in subsection 3.2. Also in this 
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Figure 4: T7ie ratto 3 A^ D = R^ ^(t' , t), for L s = 7, as a function of t' , at t = 8. 



case, in order to compute the matrix element of the chromo-magnetic operator, we 
used either R ff g(t',t) or R ff .g(t) = Ya'=o *)• The UKQCD collaboration [§ 

also used the quantity R ff g{t). In this case we make a linear fit of the form 

i^. d (t) = a 2 + 3A^ are i. (30) 



We now present our results. In fig. || we show 3\2 D = Rg_g(t',t) as a function 
of i', at fixed t = 8, obtained using the double-smeared sources Jg(x), for Ls = 7, 
at the value of the light-quark Wilson parameter K = 0.1432. In fig. [| we also give 
Rg g(t) as a function of t, for the same sources, at the same value of the light-quark 
mass. Notice that the large statistics has allowed us to perform the fit up to rather 
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Figure 5: The ratio a 2 + 3 X^H = R^,g(t), for L s = 7, as a function of t. 

large values of t (as large as t = 12), thus reducing possible systematic uncertainties 
coming from higher excitations. For comparison the UKQCD collaboration fit g(t) 
in the interval 2 < t < 5. As in the case of Aj" 116 , we use the result obtained for a 
smearing size Ls = 7 with double-smeared sources to determine the central value of 
A^ 8,16 , and take into account the variations of the result with the smearing size and 
time intervals to estimate the error. The main results are given in table |j. The value 
of A!> are is consistent, within the statistical errors, with A^ 6 being independent of the 
light-quark mass. 

In the case of A2 are we observe a slightly larger dependence of the results on the 
smearing size. We estimate the error due to the imperfect isolation of the lightest state 
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3A barc from R^yj) and Rgg{t) for L s = 7 


K 


t = 8, if = 3-5 


Fit with 6 < t < 12 


0.1425 


0.0297(16) 


0.031(3) 


0.1432 


0.0297(18) 


0.031(4) 


0.1440 


0.030(2) 


0.030(4) 


0.145431 


0.030(3) 


0.029(5) 



Table 4: Values of 3 e obtained as explained in the text. 

to be 0.005. Thus we quote 

3 A!f rc = 0.030 ±0.003 ±0.005. (31) 

In order to obtain the matrix element of the renormalized operator, we have evaluated 
the constant Z- a using the one-loop perturbative calculation of ref. || and the same 
boosted coupling as in subsection |3.3j . This gives 

Z. d = 1.85, (32) 

a very large correction that suggests that higher-order contributions can be important 
(even with the standard definition of the lattice strong coupling constant the one-loop 
correction, corresponding to a value of Z- q = 1.51, is rather large). 
Using the same calibration of the lattice spacing as before, we get 

A 2 = 0.070 ± 0.015 GeV 2 , (33) 

where the error combines in quadrature the statistical error, the systematic errors dis- 
cussed before, and the error from the calibration of the lattice spacing. Our result, 
which corresponds to a mass splitting M|, -M| = A\ 2 = 0.280±0.060 GeV 2 , although 
in agreement with the calculation of refs. pi |33|, is about a factor of 2 smaller than 
the experimental result. One may wonder whether the discrepancy originates from the 
fact that the one-loop perturbative corrections needed to obtain the matrix element of 
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the renormalized operator are so large [], or whether it is due to the quenched approx- 
imation, or to a physical reason. Notice, however, that QCD sum rule determinations 
of the splitting have given results that are in agreement with the experimental value 
[21, [2^]. We believe, however, that the possibility that our result for A2 is correct, and 
that higher-order corrections in 1/mg are still important for the B*-B mass difference, 
remains open. Using Wilson fermions for propagating quarks, one also obtains hyper- 
fine splittings which are smaller than the experimental numbers. This is a different 
problem, however, which is related to the presence of a spurious chromo-magnetic term 
of 0(a) present in the lattice action. 

From table |4| we observe a very mild dependence of A2 on the mass of the light 
quark, which would make M|, — Mg slightly larger than Mj|* — M|. 



5 Conclusions 

We have presented the results of a high-statistics lattice calculation of the kinetic 
energy Ai and of the matrix element A2 of the chromo-magnetic operator of the heavy 
quark in a 5-meson in the HQET. 

The results for Ai significantly improve a previous calculation presented in ref. M. 
Our best estimate is Ai = 0.09 ±0.14 GeV 2 . In order to make a meaningful comparison 
of this result with other theoretical determinations of the same quantity, a perturbative 
calculation of the relation between the different definitions of the renormalized kinetic 
energy operator is needed, at a sufficient degree of accuracy JTOj] . This calculation is 
missing to date, since both on the lattice and in the continuum only one-loop results 
are known. 

For A2 we obtain the value 0.070 ± 0.015 GeV 2 , which corresponds to Mg, — Mg = 
0.280 ± 0.060 GeV 2 , which is about half of the experimental number M|» - iW| = 
0.485 GeV 2 . This discrepancy, which is common to all lattice results, could have several 
8 A non-perturbative calculation of Z g q using the method proposed in ref. j?J would help to clarify the 
situation. 
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reasons, e.g. the renormalization of the relevant operator, the quenched approximation, 
or higher-order terms in the heavy-quark expansion, and needs to be clarified. 
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